Resonance Superfluidity: Renormalization of Resonance Scattering Theory 
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We derive a theory of superfluidity for a dilute Fermi gas that is valid when scattering resonances 
are present. The treatment of a resonance in many-body atomic physics requires a novel mean- 
field approach starting from an unconventional microscopic Hamiltonian. The mean-field equations 
incorporate the microscopic scattering physics, and the solutions to these equations reproduce the 
energy-dependent scattering properties. This theory describes the high-T c behavior of the system, 
and predicts a value of T c which is a significant fraction of the Fermi temperature. It is shown that 
this novel mean-field approach does not break down for typical experimental circumstances, even 
at detunings close to resonance. As an example of the application of our theory we investigate the 
feasibility for achieving superfluidity in an ultracold gas of fermionic 6 Li. 
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I. INTRODUCTION 

The remarkable accomplishment of reaching the regime 
of quantum degeneracy [lj] in a variety of ultracold atomic 
gases enabled the examination of supcrfiuid phenomena 
in a diverse range of novel quantum systems. Already 
many elementary aspects of superfluid phenomena have 
been observed in bosonic systems including vortices 0. 
The challenge of achieving superfluidity in a Fermi gas re- 
mains, however, although it appears possible that this sit- 
uation may change in the near future. A number of can- 
didate systems for realizing superfluidity in a fermionic 
gas appear very promising and it is currently the goal 
of several experimental efforts to get into the required 
regime to observe the superfluid phase transition. So far 
both fermionic potassium |^] and lithium [^J^] have been 
cooled to the microkelvin regime and are well below the 
Fermi temperature by now — a precursor step for super- 
fluidity. 

In order to make the superfluid phase transition exper- 
imentally accessible, it will likely be necessary to utilize 
the rich internal hyperfine structure of atomic collisions. 
Scattering resonances, in particular, may prove to be ex- 
tremely important since they potentially allow a signifi- 
cant enhancement of the strength of the atomic interac- 
tions. It is anticipated that by utilizing such a scattering 
resonance one may dramatically increase the critical tem- 
perature at which the system becomes unstable towards 
the formation of Cooper pairs, thus bringing the critical 
temperature into the experimentally accessible regime. 

In spite of its promise, this situation poses a number 
of fundamental theoretical problems which must be ad- 
dressed in order to provide an adequate minimal descrip- 
tion of the critical behavior. The scope of the complex- 
ities that arise in treating a scattering resonance can be 
seen by examining the convergence of the quantum ki- 
netic perturbation theory of the dilute gas. In this theory 
the small pa rame ter is known as the gaseous parameter; 
defined as V 



where n is the particle density and a 



is the scattering length. Formally, when the scattering 
length is increased to the value at which na 3 ~ 1, conven- 
tional perturbation theory breaks down §J^]. This situ- 
ation is commonly associated with the theoretical treat- 
ment of strongly interacting fermionic systems where 
high-order correlations must be treated explicitly. 

In this paper, we show that an unconventional mean- 
field theory can still be appropriately exploited under 
the condition that the characteristic range R of the po- 
tential is such that nR 3 <C 1 (while na 3 > 1). The core 
issue is that around a resonance, the cross-section be- 
comes strongly dependent on the scattering energy. This 
occurs when either a bound state lies just below thresh- 
old, or when a quasi-bound state lies just above the edge 
of the collision continuum. In both cases, the scattering 
length — evaluated by considering the zero energy limit 
of the scattering phase shift — does not characterize the 
full scattering physics over the complete energy range of 
interest, even when in practice this may cover a range of 
only a few microkelvin. 

The paper is outlined as follows. In Section II, we 
present a systematic derivation of the renormalized po- 
tentials for an effective many-body Hamiltonian. This 
requires a detailed analysis of coupled-channels scatter- 
ing. In Section III, we derive the resonance mean-field 
theory. In Section IV, we present the thermodynamic 
solutions allowing for resonance superfluidity. We apply 
our theory to the specific case of 6 Li and determine the 
critical temperature for the superfluid phase transition. 
In Section V, we consider the validity of the mean field 
approach in the case of resonance coupling, and establish 
the equivalence with previous diagrammatic calculations 
of the crossover regime between fermionic and bosonic 
superconductivity. 
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II. TWO-BODY RESONANCE SCATTERING 

The position of the last bound state in the interatomic 
interaction potentials generally has a crucial effect on 
the scattering properties. In a single-channel system, 
the scattering process becomes resonant when a bound 
state is close to threshold. In a multi-channel system the 
incoming channel (which is always open) may be cou- 
pled during the collision to other open or closed channels 
corresponding to different spin configurations. When a 
bound state in a closed channel lies near the zero of the 
collision energy continuum, a Feshbach resonance |l may 
occur, giving rise to scattering properties which are tun- 
able by an external magnetic field. The tuning depen- 
dence arises from the magnetic moment difference Afi mas 
between the open and closed channels |9| . This gives rise 
to a characteristic dispersive behavior of the s-wave scat- 
tering length at fields close to resonance given by 



1 - 



AB 

B — Br 



(1) 



where Ob g is the background value which may itself de- 
pend weakly on magnetic field. The field-width of the 
resonance is given by AB, and the bound state crosses 
threshold at a field- value Bo. The field-detuning can 
be converted into an energy-detuning v by the relation 
v = (B — £?o)Au mag . An example of such a resonance 
is given in Fig. ol, where a coupled channels calculation 
is shown of the scattering length of 6 Li for collisions be- 
tween atoms in the (f,m f ) = (1/2,-1/2) and (1/2,1/2) 
state jl(J. The background scattering length changes 
slowly as a function of magnetic field due to a field- 
dependent mixing of a second resonance which comes 
from the triplet potential. This full coupled channels 
calculation includes the state-of-the-art interatomic po- 
tentials pd[ | and the complete internal hyperfine struc- 
ture [Bl. 
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The scattering length is often used in many-body 
theory to describe interactions in the s-wave regime. 
That the scattering length completely encapsulates the 
collision physics over relevant energy scales is implic- 
itly assumed in the derivation of the conventional 
Bardeen-Cooper-Schrieffer (BCS) theory for degenerate 
gases [ pT|JT5| , as well as the Gross-Pitaevskii description 
of Bose-Einstein condensates. However, the scattering 
length is only a useful concept in the energy regime where 
the s-wave scattering phase shift depends linearly on 
the wavenumber fc, i.e. Sq = —ka. For a Feshbach reso- 
nance system at a finite temperature there will always be 
a magnetic field value where this approximation breaks 
down and the scattering properties become strongly en- 
ergy dependent. In close proximity to a resonance, the 
scattering process then has to be treated by means of the 
energy dependent T-matrix. 

Only the exact interatomic interaction will reproduce 
the full T-matrix over all energy scales. However, since 
only collision energies in the ultracold regime (of order 
microkelvin) are relevant, a much simpler description is 
possible. If the scattering length does not completely 
characterize the low energy scattering behavior in the 
presence of a resonance, what is the minimal set of pa- 
rameters which do? 

As illustrated in Fig. ||[ we proceed to systematically 
resolve this question by the following steps. We start 
from a numerical solution of the complete coupled chan- 
nels scatte ring problem for a given real physical system. 
In Section II A we demonstrate that the results of these 



full numerical calculations can be adequately replicated 
by giving an analytic description of resonance scattering 
provided by Feshbach's resonance theory. The point of 
this connection is to demonstrate that only a few param- 
eters are necessary to account for all the collision proper- 
ties. This implies that the scattering model is not unique. 
There are many microscopic models which could be de- 
scribed by the same Feshbach theory. In Section [IB 



we show this explicitly by presenting a simple double 
well model for which analytic solutions are accessible. 
Thereby we derive a limiting model in which the range of 
the square well potentials and coupling matrix elements 
are taken to zero. This leads in Section II C to a scat- 



tering model of contact potentials. We show that such a 
scattering solution is able to reproduce well the results 
of the intricate full numerical model we began with. The 
utility of this result is that, as will be apparent later, it 
greatly simplifies the many-body theoretic description. 
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FIG. 1. Scattering length as a function of magnetic field, 
for the (f,m f ) = (1/2, -1/2) and (1/2, 1/2) mixed spin chan- 
nel of 6 Li. 
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Full CC scattering 



Feshbach scattering 



Square-well scattering 



Contact scattering 



FIG. 2. Sequence of theoretical steps involved in formu- 
lating a renormalized scattering model of resonance physics 
for low energy scattering. The starting point is a full coupled 
channels (CC) calculation, which leads via an equivalent Fes- 
hbach theory, and an analytic coupled square well theory, to 
a contact potential scattering theory which gives the renor- 
malized equations for the resonance system. 



A. Feshbach resonance theory 

Here we briefly describe the Feshbach resonance for- 
malism and derive the elastic S- and T-matrices for two- 
body scattering. These matrices arc related to the tran- 
sition probabilities for scattering from an initial channel 
a to a final channel [3. A more detailed treatment of this 
formalism can be found in the literature ||. 

In Feshbach resonance theory two projection opera- 
tors P and Q are introduced which project onto the sub- 
spaces V and Q. These subspaces are two orthogonal 
components which together span the full Hilbert space of 
both scattering and bound wavefunctions. The open and 
closed channels are contained in V and Q, respectively. 
The operators P and Q split the Schrodinger equation 
for the two-body problem into two parts: 



(E-H PP )\ij p ) =H PQ \^) 
(E-H QQ )\^) =H QP \i> p ) 



(2) 
(3) 



where Hpp — PHP, Hpq = PHQ, etc., and ip is the 
total scattering wavefunction. The projections on the 
two sub-spaces are indicated by P\tjj) = |VO an( l Q\fp) = 
\ip < 3). The Hamiltonian H = Ho + V consists of the sum 
of the single-particle interactions Hq and the two-body 
interaction V. Eq. (||) can be formally solved 



|V»°) = 



H 



QQ 



(4) 



where E + = E + i5 with 6 approaching zero from positive 
values. Substituting this result into Eq. (||), the open 
channels equation can be written as (E — H c s)\tp p ) = 0, 
where 

ffeff = Hpp + Hpq— — H QP . (5) 

- tlQQ 

The resolvant operator is now expanded in the discrete 
and continuum eigenstates of Hqq : 

HpQ\(t>i){(j>i\HQp 



Hqs = H 



pp 



H PQ \<j>{e)){cj>{e)\H QP 
E+ -e 



(6) 



de. 



Here the ej's are the uncoupled bound-state eigenvalues. 
In practice, only a few bound states will significantly af- 
fect the open-channel properties. In this paper, we will 
consider either one or two bound states and neglect the 
continuum expansion in Eq. (^). Then the formal solu- 
tion for \tp } is given by 

Hpq^^Hqp^) 



1 



Hpp 



£■ 



E 



(7) 

where \ip„) is the eigenstate of the direct interaction 
Hpp that satisfies the outgoing wave boundary condi- 
tion in channel a. By multiplying from the left with 
(X/?l^) where \xp) is an unscattered state in the outgo- 
ing channel (3, the left-hand-side becomes the T-matrix 
for the total scattering process. The unscattered state 
is related to the scattering wave- function |V , ( a _ ) with in- 
coming boundary conditions via 

V 



\X(3) 



\Xp) 



E — Hp p 

The T-matrix giving the transition amplitude is then 



(8) 



T — T p 
1 8a — 1 /3a 



£ 



(^-\H pq \^)(^\H Q p\^ p ) 
E-e, 



(9) 



where T p a is the amplitude for the direct (non-resonant) 
process. From the T-matrix we can easily go to the 
S'-matrix that is defined as Sp a = (ipg\ipa)- Since 
we consider s-wave scattering only, in our case there 
exists a simple relation between the S- and T-matrix: 
Sf3 a = 1 — 27ri7g a p6| , and this allows us to rewrite 
Eq. (§) as 



S/3a — Sp a 



_ „ 2-Ki{i ) + \H PQ \<t> i ){<t> i \H QP \^ 1 
2^ b P-y 2^ 

7 i 



E - 



(10) 

The non-resonant factors Sp describe the direct scat- 
tering process from an open channel 7 to the outgoing 
channel j3. Returning to Eq. (Q), we can solve for the 
component ((f>i\HQ P \tjj p ) by multiplying both sides with 
(4>i\Hnp- 
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1. Single resonance 

For the case of only one resonant bound state and only 
one open channel, the solution of Eq. (Q) gives rise to the 
following elastic S-matrix element (we will omit now the 
incoming channel label a): 



S = S' 



1 - 



2m\{^ p +\H PQ \ 



E-e x - (<h\HQP E±± Hpp BpQ\<h) 



(11) 



The non-resonant S'-matrix is related to the background 
scattering length via S p — exp[— 2ikdb s ]. The term in 
the numerator gives rise to the energy-width of the res- 
onance, r = 2 / K\(tp p+ \HpQ\(j>i)\ 2 7 which is proportional 
to the incoming wavenumber k and coupling constant g~\ 
|l7f . The bracket in the denominator gives rise to a shift 
of the bound-state energy, and to an additional width 
term iT/2. When we denote the energy-shift between 
the collision continuum and the bound state by v\, and 
represent the kinetic energy simply by H 2 k 2 /m, the S'- 
matrix element can be rewritten as 



S(k) = 



—Uka^rj 



2ik\g x \ 



^ 1 -^)+ik\- 9l \ 



(12) 



The resulting total scattering length has exactly the dis- 
persive lineshape for the resonant scattering length which 
we presented originally as Eq. ([[]). 



2. Double resonance 

Often more than one resonance may need to be con- 
sidered. For example, the scattering properties for the 
(1/2,-1/2) + (1/2,1/2) channel of 6 Li are dominated 
by a combination of two resonances: a triplet poten- 
tial resonance and a Feshbach resonance. This can be 
clearly seen from Fig. [I], where the residual scattering 
length, which would arise in the absence of the Fesh- 
bach resonance coupling, would be very large and nega- 
tive and vary with magnetic field. This can be compared 
with the value of the non-resonant background scatter- 
ing length for the triplet potential for Li which is only 
31 ao, which is an accurate measure of the characteristic 
range of this potential. An adequate scattering model for 
this system therefore requires inclusion of both bound- 
state resonances. Since for 6 Li the coupling between 
these two bound states is small, it will be neglected in 
the double resonance model presented here. The double- 
resonance S-matrix, with again only one open channel, 
follows then from Eq. ( |l0| ) and includes a summation over 
two bound states. After solving for the two components 
(<f>i\HQp\ip p ) of wave function \ip ), the S-matrix can be 
written as 



S(k) 



2ik(\g 1 \ 2 A 2 + \g 2 \ 2 A 1 ) 



ik(\g 1 \ 2 A 2 + \g 2 \ 2 A 1 )-A 1 A 2 



(13) 



with Ai = (p\ — h 2 k 2 /m)4:TTh 2 /m, where v\ and g\ are 
the detuning and coupling strengths for state 1. Equiva- 
lent definitions are used for state 2. Later we will show 
that this simple analytic Feshbach scattering model mim- 
ics the coupled channels calculation of 6 Li. The parame- 
ters of this model, which are related to the positions and 
widths of the last bound states, can be directly found 
from a plot of the scattering length versus magnetic field 
as given, for example, by Fig. El. The scattering length 
behavior should be reproduced by the analytic expression 
for the scattering length following from Eq. (13): 



a bg 



m 
inh 2 



(14) 



The advantage of a double-pole over a single-pole S- 
matrix parametrisation is that we can account for the 
interplay between a potential resonance and a Feshbach 
resonance, which in principle can radically change the 
scattering properties. This interplay is not only im- 
portant for the description of 6 Li interactions, but also 
for other atomic systems which have an almost resonant 
triplet potential, such as bosonic 133 Cs [[18 19 



and 85 Rb 



In the many-body part of this paper, Section III, the 



scattering properties are represented by a T-matrix in- 
stead of an S-matrix. We have shown in the above that 
in our case there exists a simple relation between the two, 
however, the definition for T in the many-body theory 
will be slightly different in order to give it the conven- 
tional dimensions of energy per unit density: 



T(k) 



2nh 2 i 
mk 



[S(k) - 1] 



(15) 
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B. Coupled square- well scattering 
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FIG. 3. Illustration of the coupled square well sys- 
tem. Outer region r > R: the solid line corresponds 
to the open channel potential P, and the dotted line to 
the closed channel potential Q. The wave functions are 
given by up(r) ~ sin kpr and uq(t) ~ exp(— kQr), respec- 
tively. Inner region r < R: the solid and dotted lines 
correspond to the molecular potentials Vi and V2, respec- 
tively. The wave functions are given by lii(r) ~ sinfcir and 
U2{r) ~ sin(A;2r). The dashed line corresponds to the kinetic 
energy E in the open cha nnel. The wavevectors are defined as: 
k P = VmE/H, kg = ^m(e - E)/h, ki = \/m(E + Vi)/h, 

k2 = v/ + V2 — e)/h. The detuning e can be chosen such 
that a bound state of square-well potential Vi enters the col- 
lision continuum, causing a Feshbach resonance in the open 
channel. 

In this subsection we describe the coupled-channels ex- 
tension of a textbook single-channel square-well scatter- 
ing problem. One reason that this model is interesting 
to study is because we can take the limit of the potential 
range R — > 0, thus giving an explicit representation of a 
set of coupled delta-function potentials which simplifies 
the description in the many-body problem to follow. 

The scattering equations for such a coupled system are 
written as 



center of mass frame. The coupled square well model en- 
capsulates the general properties of two-body alkali in- 
teractions. There we can divide the internuclear sepa- 
ration into two regions: the inner region where the ex- 
change interaction (the difference between the singlet and 
triplet potentials) is much larger than the hyperfinc split- 
ting, and the outer region where the hyperfine interaction 
dominates. Here we make a similar distinction for the 
coupled square wells. In analogy to the real singlet and 
triplet potentials, we use for the inner region two artificial 
square- well potentials labeled as V\ and V 2 . We take the 
coupling <?(r) to be constant over the range of the square- 
well potentials r < R, and to be zero outside this range 
(see Fig. [|). Then the problem can be simply solved by 
means of basis rotations at the boundary R giving rise 
to simple analytic expressions. For r > R, we therefore 
consider one open channel and one closed channel, with 
wavenumbers kp and kQ . In analogy with a real physical 
system, we can refer to the inner range channels (r < R) 
as a molecular basis, and the channel wave functions are 
just linear combinations of the U\ and u 2 wave functions. 
At the boundary i?, these wave functions have accumu- 
lated a phase <fii — k\R and <f> 2 = k 2 R. The coupling 
strength is effectively given by the basis-rotation angle 9 
for the scattering wave functions: 



u P (R) 
u Q {R) 



cos a — sin ( 
sin 9 cos 9 



ui{R) 
u 2 (R) 



(18) 



allowing for an analytic solution of the scattering model. 
This leads to the following expression for the S'-matrix: 



S=e~ MKpU [l- (19) 
(—2ikp{k 2 cot 4>2 cos 2 6 + kQ + ki cot 4>\ sin 2 8))/ 
(kpkQ + ki cot 4>i (kp sin 2 9 — kQ cos 2 9) + 



ik 2 cot 4>2{k\ cot 4>\ + kp cos 2 9 + kQ 



'O))] 



EiP Q (v) = 



__v +y» 

m 



V Q (r) 



^ p (r)+ 5 (r)^ Q (r), 

^(r)+ 3 *(r)V> P (r), 



with e the energy-shift of the closed channel with respect 
to the collision continuum and E = h 2 k 2 /m the rela- 
tive kinetic energy of the two colliding particles in the 



An extention to treat more than two coupled potentials, 
which would be required to model more than one reso- 
nance, is also straightforward. 

The parameters of the two wells have to be chosen 
such that the results of a real scattering calculation are 
reproduced for a given physical system. In fact all the 
parameters are completely determined from the field de- 
pendence of the scattering length, and all other scattering 
properties, such as the energy-dependence of the scatter- 
ing phase shift, can then be derived. First we choose a 
range R, typically of the order of an interatomic potential 
range (100 do) or less. Now we have only to determine 
the set of parameters V\, V 2 and 9. The potential depth 
(16) Vi is chosen such that the scattering length is equal to the 
background scattering length ab g , while keeping 9 = 0. 
Also, Vi should be large enough that the wavenumber 
^ 'k\ depends weakly on the scattering energy. Then, we 
set 9 to be non-zero, and change the detuning until a 
bound state crosses threshold, giving rise to a Feshbach 
resonance. The value of V 2 is more or less arbitrary, but 
we typically choose it to be larger than V\. Finally, we 
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change the value of 8 to give the Feshbach resonance the 
desired width. 

We will later show that the resulting scattering prop- 
erties converge for R — ► 0. In Fig. [| the coupled square- 
well system is compared with the Feshbach scattering 
theory, for 40 K scattering parameters. Even despite the 
fact there is a strong energy-dependence of the T-matrix, 
the two scattering representations agree very well. 




EOiK) 



FIG. 4. Comparison of the real part of the T-matrix for 
coupled square-well scattering (solid line) with a potential 
range R = lao, to Feshbach scattering (dashed line), for a 
detuning that yields a scattering length of about -2750 ao- 
Similar good agreement is found for all detunings. 



C. Contact potential scattering and renormalization 

In this subsection the Lippmann-Schwinger scattering 
equation is solved for a resonance system with contact 
potentials. As in the previous subsection, we make use 
of an open subspace that is coupled to a closed subspace. 
The contact potentials are defined by: 



V p (r) = V p 5(r), 
V Q (r) = V Q S(r), 
g(r) =g5(r), 



(20) 
(21) 



with S(r) is the three-dimensional Dirac delta-function. 
Here V p (r) is the open channel potential with strength 
V p . The function V^r) is a closed channel potential 
with strength V®, and g(r) is a coupling between the 
closed and open channel with strength g. The procedure 
of renormalizatio n rel ates the physical units (abg, <?i, and 



Vi) from Section |II A| to these parameters of the contact 
potential scattering model for given momentum cutoff; 
a relationship for which we will now obtain explicit ex- 
pressions. The first step is to solve again the scattering 
Eqns. ( |l6|) and (|l7| ) for these contact potentials. As we 
have seen in subsection [I A , we can formally solve the 
bound-state equations, and make use of Eq. (^3j) to ex- 
pand the Green's function in bound-state solutions. In 
this case it can be written as 



<A Q (r)=E 



^(r)/dV0?V)5*(r')^V) 



E 



(22) 



with <j)~f (r) a bound state solution and e, its eigenenergy. 
We now define an amplitude for the system to be in this 
bound state that will later be useful in the mean-field 
equations: <j>i = ((f)f\ip®), and together with the open 
channel equation and the definition <7i(r) = g(r)cf>f(r), 
we get a new set of scattering equations 



2l2 



h 4 k 



-—V 2 r + V p (r) 



V/>) + 5>(r)^, (23) 



2 1.2 



(24) 



The energy-difference between the bound-state energy 
and the threshold of the collision continuum is given by 
Vi. The open channel solution for Eq. (^3|) can be formu- 
lated as 



V> P (r) = x(r) 



Ak\r—r' I 



4irh 2 



d ^ uT-JT [V P (r>)^ P (r>) 



(25) 



X« +/(*)- 



. jkr 



CO. 



Here x( r ) is the unscattered wavefunction, and in the 
other term we recognize the scattered part that is usu- 
ally formulated in terms of the scattering amplitude f(6). 
The momentum representation of this last line is M: 



V- P (p) = (27r) 3 J(k-p) 



iir f(jk,p) 



k 2 



P 



id 



(26) 



Combining Eq. (|26|) with our expression for the scattering 
amplitude we find 



47rr 



f(k,k') = V p + 



1 



(27T) 



rV 1 



Airh z 



d 3 P 



f(k,P) 



h 2 P 2 



+ i5 
(27) 



The typical temperature range of a system we are in- 
terested in will only allow for elastic s-wave scattering, 
therefore the scattering amplitude has no angular de- 
pendence, and incoming and outgoing wavenumbers are 
the same, i.e. k = k' . The scattering amplitude can 
then be simply linked to the T-matrix via the relation 
T(k) = — (inh 2 1 to) / (fc) . The integral has a principal- 
value part, and the integration ranges from zero to a 
momentum cutoff K. Eq. (B7j) then has as solution 



T{k) = V^ - 



p VPm T(k) 



2n 2 h 



2*2 



, k vk , 

K — arctanh 1 k 

K 2 



(28) 



G 



This is a variant of the Lippmann-Schwinger equation. 
The closed channel scattering solutions are now used to 
eliminate the amplitude functions fa. In Fourier space, 
Eq. has the form 



2l2 



h z k 



-fa = Vifa + g* 



1 



(2tt) 3 



(29) 



After substitution of Eq. (|26|) the expression for 
linked to the T-matrix: 



is 



1 



T(k) [K - arctanh 



*]) 



h 2 k 2 



(30) 



Eliminating fa from Eq. (|2§|) gives a complete expression 
for the Lippmann-Schwinger equation 



T(jfc) = V 



V 171 



T(k) 



k in 
K — arctanh— + — k 
K 2 



2tt 2 ^ 2 

M 1 - 2^f T ( fc ) l K ~ arctanh^ 



:*]) 
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(31) 



Similar to the Feshbach and coupled square- well prob- 
lems, the k — > behavior of T(fc) should reproduce the 
scattering length, and, the result should not depend on 
the arbitrary momentum cutoff K. For an analytic ex- 
pression of the scattering length, we conveniently use the 
Feshbach representation. A comparison between the lat- 
ter and the expression for the scattering length a that 
results from solving Eq. (|3l]), tells us how to relate the 
coupling constants for contact-scattering to the Feshbach 
coupling constants. By making use of the definitions 
r = (l-aU)' 1 , a = mK/(2n 2 Ti), and U = Anh 2 : a hg /m, 
we find the very concise relations 



V 1 



ru, 



(32) 



which is valid also in the case where no resonance is 
present, and in addition 



9i 



vi + agi9i- 



(33) 
(34) 



for the open channel potential and the first resonance. 
For the second resonance, if present, we find 



.92 



f)2 



agl/vx+T- 1 ' 
v 2 = V2 + ag 2 g2- 



(35) 
(36) 



Obviously, our approach an be systematically extended 
further, order by order, to give an arbitrarily accurate 
representation of the microscopic scattering physics. 

These expressions we refer to as the renormalizing 
equations of the resonance theory since they remove the 
ultra-violet divergence which would otherwise appear in 
the field-equations. Any many-body theory based on con- 
tact scattering around a Feshbach resonance will need to 



apply these expressions in order to renormalize the the- 
ory. These equations (^)-(^) therefore represent one of 
the major results of this paper. 

In Fig. ^| the T-matrix as a function of energy is shown 
for contact scattering, in comparison with the square- well 
scattering for different values of the potential range. The 
contact-scattering model is demonstrated to be the limit- 
ing case of the coupled square well system when R — > 0. 
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FIG. 5. Comparison of the real part of the T-matrix for 
coupled square- well scattering for three different values of the 
potential range: R = 100ao (dash-dotted line), R = 30ao 
(dashed line), R = lao (solid line). The interaction param- 
eters for 40 K have been used here, and the magnetic field 
is chosen such that a scattering length of a — 300ao is ob- 
tained. Also plotted is the T-matrix for contact scattering, 
which clearly agrees very well as it coincides with the solid 
line of the double-well scattering. 



D. Discussion of different models 



In subsection II C it has been shown that the resonance 



contact scattering representation is the limiting case of 
the coupled square-well system, when the range of t he 
potentials is taken to zero. Also, in subsection II B it 
has been shown that the double-well system is in good 
agreement with the Feshbach scattering theory. Now we 
will show how well these scattering representations agree 
with the full numerical coupled channels calculation p0[ . 
In Fig. |^ we show the real and imaginary parts of the T- 
matrix applied to the case of 6 Li, and compare the cutoff 
and Feshbach scattering representations to a full coupled 
channels calculation. The agreement is surprisingly good, 
and holds basically for all magnetic fields (i.e. similar 
agreement is found at all detunings). 
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FIG. 6. (a) Real part of the T-matrix as a function of col- 
lision energy, for the Feshbach model and the cutoff model 
(overlapping solid lines), and for a coupled channels calcula- 
tion (dashed line). The atomic species considered is 6 Li, for 
atoms colliding in the (1/2,-1/2) + (1/2, 1/2) channel, (b) 
Same as (a) for the imaginary part. 

In this section we have discovered a remarkable fact 
that even a complex system including internal structure 
and resonances can be simply described with contact po- 
tentials and a few coupling parameters. This was known 
for off resonance scattering where only a single parameter 
(the scattering length) is required to encapsulate the col- 
lision physics at very low temperature. However, to our 
knowledge this has not been pointed out before for the 
resonance system, where an analogous parameter set is 
required to describe a system where the scattering length 
may even pass through infinity. We have shown in a very 
concise set of formulas how to derive the resonance pa- 
rameters associated with contact potentials. This result 
is important for the incorporation of the two-body scat- 
tering in a many-body system, as we will show later in 
this paper. 

Other papers have also proposed a simple scattering 
model to reproduce coupled channels calculations |^l|,^2| . 
In these papers real potentials are used, and they give 



a fair agreement. Here, however, we use models that 
need input from a coupled channels calculation to give 
information about the positions of the bound states and 
the coupling to the closed channels. All this information 
can be extracted from a plot of the scattering length as 
a function of magnetic field. 



III. MANY-BODY RESONANCE SCATTERING 

We will now proceed to a many-body description of res- 
onance superfluidity and connect it to our theory of the 
two-body scattering problem described earlier. This sec- 
tion explains in detail the similar approach in our papers 
devoted to resonance superfluidity in potassium p3| , p4j . 
The general methods of non-equilibrium dynamics has 
been described in |25| and we have applie d them in the 
context of condensed bosonic fields |26|j27f| . 

In the language of second quantization, we describe 
the many-body system with fermionic fields ?p a (x) which 
remove a single fermionic particle from position x in in- 
ternal electronic state a, and molecular bosonic fields 
4> i (x) which annihilate a composite bound two-particle 
excitation from space-point x in internal configuration i. 
These field operators and their adjoints satisfy the usual 
fermionic anti-commutation rules 

{^(^O.V^C 3 ^)} = 5{xi - X 2 )S ai<T2 = 6l2, 

and bosonic commutation rules 



o. 



(38) 



respectively. In here and the following discussion, we will 
also try to simplify the notational complexity by adopting 
the notation convention of many-particle physics. This 
means, we will identify the complete set of quantum num- 
bers uniquely by its subscript index, i.e., {x\, <t\} = 1. If 
only the position coordinate is involved, we will use bold 
face X2 = 2. 

In the double resonance case of lithium, we have to 
distinguish only two internal atomic configurations for 
the free fermionic single particle states a = {T-l} an d 
we need at most two indices i — {1,2} to differentiate 
between the bosonic molecular resonances. 

The dynamics of the multi-component gas is governed 
by a total system Hamiltonian H = Hq + H±, which 
consists of the free evolution Hamiltonian Hq and the in- 
teractions H\ between atoms and molecules. We assume 
that the free dynamics of the atoms and molecules is de- 
termined by their kinetic and potential energies in the 
presence of external traps, which is measured relative to 
the energy /i of a co-rotating reference system. Thus, we 
define 
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H a (x) 



2m 



v 2 + y CT (x)-/i, 



(39) 
(40) 



Here, m denotes the atomic mass as used previously, 
M = 2 to is the molecular mass, /i TO = 2 fj, is the energy 
offset of the molecules with respect to the reference sys- 
tem, V a (x) are external spin-dependent atomic trapping 
potentials, and V l m (x) are the external molecular trap- 
ping potentials. The molecular single-particle energy has 
an additional energy term Vi that accounts for the de- 
tuning of the molecular state i relative to the threshold 
of the collision continuum. 

The binary interaction potential V p {x\ — X2) accounts 
for the non-resonant interaction of spin-up and spin-down 
fcrmions, and coupling potentials gi{x\ — x-i) convert free 
fermionic particles into bound bosonic molecular excita- 
tions. Thus, we find for the total system Hamiltonian of 
the atomic and molecular fields: 



H — Ho + Hi, 



(41) 



where the free Hq and interaction contributions Hi are 
defined as 



dl^ 7/^(1)^(1)^(1) 
a 

dlJ2 #(1)^(1)^(1), (42) 

i 

dld2 {$(1)$(2)V P (1 - 2)^(2)7^(1) 

)g*(l- 2)^(2)^(1) +H.c. J. 

(43) 



H a = 
+ 

Hi = 

+E 



Here, H.c. denotes the Hermitian conjugate. In the 
present picture, we deliberately neglect the interactions 
among the molecules. Several other papers have treated 
a Feshbach resonance in a related way p8|-|3"H . 

In order to derive dynamical Hartree-Fock-Bogoliubov 
(HFB) equations from this Hamiltonian, we also need to 
define a generalized density matrix to describe the state 
of the fermionic system |3^] and an expectation value for 
the bosonic molecular field. The elements of the 4x4 
density matrix Q are given by 

g m {12) = {A\{x 2 )A p {xi)), (44) 

A(x)=(j)^{x),$ i (x),^x),$\{x)} , (45) 

and symmetry broken molecular fields are defined as 

i (l) = (^(x 1 )). (46) 

As usual, we define the quantum averages of an arbitrary 
operator (?) with respect to a many-body density matrix p 



by (<Z)) = Tr[(?)p], and we calculate higher order correla- 
tion functions by a Gaussian factorization approximation 
known as Wick's theorem |52| . The structure of the 4x4 
density matrix 



0(12) 



ftv(12) ^(12) 
-&a(12)* 5i2-0at(12) 4 



(47) 



is very simple, if one recognizes that it is formed out of a 
2x2 single particle density matrix Qjj , a pair correlation 
matrix Qa and obviously the vacuum fluctuations (5i2- 
The single particle submatrix is given by 



M12) 



0„ T (12) £ m (12) 

g m {2iy 0,4(12) 



(48) 



where Q na (12) = (^(x2)ifj a (xi)) is the density of spin- 
up and down particles and C* m (12) = (^(x2)tp^(xi)) 
denotes a cross-level coherence, or "magnetization" be- 
tween the states. The pair-correlation submatrix Qa is 
defined analogously as 



£U(i2) 



0„ t (i2) g p (i2) 
-g p {2i) g al (i2) 



(49) 



where Q acr {12) — (ip (T (x2)il> (T (xi)) is an anomalous pair- 
ing field within the same level and the usual cross- 
level pairing field of BCS theory is defined in here as 
g p (12) = (^(aB 2 )^ T (aJi)>. 



A. General dynamic Hartree-Fock-Bogoliubov 
equations of motion 

From these physical assumptions about the system's 
Hamiltonian Eq. ( pd| ) and the postulated mean fields <fii 
of Eq. ( [46"1) and g of Eq. (|47| ) , one can now derive kinetic 
equations for the expectation values (O) for an operator 
O by a systematic application of Heiscnberg's equation 



ih—O = [0,H] 



(50) 



and Wick's theorem. 

The first order kinetic equation for the Hermitian den- 
sity matrix g has the general form of a commutator and 
the time-evolution is determined by a Hermitian self- 
energy matrix £ = £° + £ . In general, one finds 

ih-g(13)= / d2 [£(12) £(23) -a(12)E(23)], (51) 
ift^(3)=flf(3)^(3) (52) 
dld2 5(±±^-3)g*(l-2)g p (12). 

First, the free evolution £° is obviously related to the 
single particle Hamiltonians of Eq. (^). In complete 
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analogy to the generalized density matrix, it has a simple 
4x4 structure 



4>i{x) = <f>i(o) = 

G(x 1 ,x 2 ) = G(xi - x 2 ) = G(r). 



(59) 
(60) 



^ ~ I -2^(12)* 



which can be factorized into 2x2 submatrices as 



£at(12) - s 12 ( H ^ 1} 



(53) 



(54) 



Secondly one obtains from the interaction Hamiltonian 
of Eq. (M3) the first order self-energy S 1 as 



E 1 (12) 



2^(12) 2^(12) 
-2^(12)* -2^(12)* 



(55) 



This assumption implies also that there can be no exter- 
nal trapping potentials present, i.e. V a {x) — V j m (x) = 0, 
as this would break the translational symmetry. 

Furthermore, we want to consider a special situation 
when there is no population difference in spin-up and 
spin-down particles G n (j" = \f\) = Gnn(r), there excists 
no cross- level coherence or "magnetization" Gm{r) = 0, 
and the anomalous pairing field Ga(r) = 0. It is impor- 
tant to note that this special scenario is consistent with 
the full evolution equation and, on the other hand, leads 
to a greatly simplified sparse density matrix: 



The normal potential matrix Ej^- has the usual structure 
of direct contributions [i.e. local Hartree potentials pro- 
portional to 612] an d exchange terms [i.e. non-local Fock 
potentials proportional to V p (l — 2)]: 



G(12) 



( Gn{r) G P (r) \ 

G n (r) -G P {r) 

-g;(r) 5(r)-Gn(r) 

\G* p (r) S(r)-Gn(r)J 

(61) 



2^(12) 



d4V p (2 - 4) 



<*120nl(44) 

Si4G m (21)* 



-<5l4<?m( 12 ) 

5i2S„ T (44) 



where r = |r| = |l — 2|. Similarly, 
tionally invariant self-energy 2(12) 



(56) 



The zeros that appear in the diagonal of the anomalous 
coupling matrix 



5m(12) 





-A(21) 



A(12) 




(57) 



reflect the fact that there is no low energy (s-wave) in- 
teraction between same spin particles due to the Pauli 
exclusion principle. The off-diagonal element defines a 
gap function as 



A(12) = V p (l - 2) G P (12) + - 2 ) <M : 



-)• 

(58) 



2(12) = 8x2 












2(1) 


-A 




-A* 


-2(1) 










-2(1) / 



one finds a transla- 
= E(l - 2) with 



(62) 



and S(a;) = — Ti /(2m)V|; — (J. + V <7 n (0) and a complex 
energy gap A = V P G P {0) + J2i 9i 0»- These assumptions 
lead to a significant simplification of the HFB equations. 

The structure of the HFB equations can be elucidated 
further by separating out the bare two-particle inter- 
actions from the many-body contributions. One can 
achieve this by splitting the self-energy into the kinetic 
energy and mean-field shifts 2 = 2° + 2 , and by sep- 
arating the density matrix into the vacuum contribution 
G° [proportional to 5(r)] and the remaining mean-fields 
G = G° + G 1 : 



B. The homogeneous limit and the contact potential 
approximation 

In this section, we will apply the general HFB equa- 
tions of motion [Eq. (|5l|)] to the case of a spatially homo- 
geneous isotropic system. Furthermore, we will approxi- 
mate the finite range interaction potentials V p (xi — x 2 ) 
and gi(xi — x 2 ) by the contact approximation as intro- 
duced in Eq. (EOh , and assume equal populations for spin- 
up and spin-down atoms. 

Spatial homogeneity implies that a physical system is 
translationally invariant. Thus, any single particle field 
must be constant in space and any two-particle quantity 
or pair-correlation function can depend on the coordinate 
difference only: 



ih-G 1 - [E ,^ 1 ] - [E\S°] = [2\S X ], (63) 



(64) 



In this fashion, we can now identify the physics of reso- 
nance scattering of two particles in vacuo [left hand side 
of Eq. flog)] from the many-body corrections due to the 
presence of a medium [right hand side of Eq. (p3|)1 

In the limit of very low densities, we can ignore many- 
bod y effe cts and rediscover Eqs. (|23| ) and ( |24j ) of subsec- 
tion II C , but given here in a time-dependent form. They 
describe the scattering problem that we have solved al- 
ready: 



0. 



(65) 
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ih-G p (r) 



+ 



m 

E 



-V"£ -2 h t + V p (r) 
9i(r)<t>i, 



G P (r) 



ift j70i = [Vi ~ Mm) <l>i + g* G P (0). 



(66) 
(67) 



The scattering solution of Eqs. ( |66[ ) and (|67() is 'sum- 
marized' by the energy-dependent two-body T-matrix, 
which we have discussed in the preceeding sections. In or- 
der to incorporate the full energy-dependence of the scat- 
tering physics, we propose to upgrade the direct energy- 
shift V p g n (r) to (T Re (k))g n (r), where (T Rc (fc)) repre- 
sents the real part of the two-body T-matrix, and (. . .) 
denotes two-particle thermal averaging over a Fermi- 
distribution. A detailed calculation of the proper up- 
grade procedure will be presented in a forthcoming pub- 
lication. 

The translationally invariant HFB Eqns. ( |63| ) and j6^ ) 
are best analyzed in momentum-space. Thus, we will 
introduce the Fourier transformed field-operators dka by 



E e 

h 



—iha 



(68) 



where is the quantization volume. If we define the 
Fourier components of the translationally invariant mean 
fields as 



g(r) = g( Xl -x 2 ) = J2 e- lk ^-^ g(k), 



(69) 



we obtain the following relations between the real-space 
density of particles n (the same for both spins) and the 
real-space density of particle pairs p 



P 



= g na (r = o) = J2e™(k) = ^Yl 

k k 

= g p (r = 0) = G P {k) = ^ E^-fci a 



fc'l / 



(70) 
(71) 



IV. THERMODYNAMICS 

In this paper we focus on the properties of thermo- 
dynamic equilibrium. Thermodynamic equilibrium can 
be reached by demanding that the grand potential $g = 
—kfrT In S at a fixed temperature has a minimal value. In 
this definition kb is Boltzmann's constant, and 5 the par- 
tition function S = Tr[exp(— -£/diag/fcf>T)]. The exponent 
containing the diagonalized Hamiltonian reads 




10 15 
Single-particle Energy (n K) 

FIG. 7. Two overlapping quasi-particle energy spectra as 
a function of single-particle kinetic energy E = h 2 k 2 /(2m). 
The total density of the gas is n = 10 14 cm -3 , the magnetic 
field is B = 900 G, and the temperature is T = 0.01T F . The 
energy-spectrum is calculated for two different values of the 
cutoff: i.e. K = 32/cf and K = The fact that the 

two lines are overlapping to the extent that the difference is 
difficult to see shows the renormalization in practice, demon- 
strating the validity of Eq. (||)-(||). 

-Hdiag = Z^jVi ~ Mm)|0i| 2 (75) 
i 

+ + E k ("L T "fe T + & H & kl - !)) » 



The Fourier-transformed HFB equations are now local in 
momentum-space 



ih^-gk) 
at 



ih- 



dt 



m),g(k)], 

(vi -y>m)4>i + g*p, 



and the self-energy is given by 



E(fc) = 











A 





s fe 


-A 








-A* 


-Sfc 





A* 








-E 



(72) 
(73) 

(74) 



In here, the upgraded single particle excitation energy 
is now = tk — A 1 + (T'k' e )n, ek = h 2 k 2 /2m denotes 
the kinetic energy and the gap energy is still defined as 
A = V p P + J2 l g l fr. 



which is a quadratic approximation to the original Hamil- 
tonian. The energy spectrum Ek results from a local di- 
agonalization by the Bogoliubov transformation of the 
self-energy matrix E^, at each fe, where the obtained 

" 2 + A 2 . Note that 



quasi-particle spectrum is Ek 



the first summation term in -ffdiag results from a contri- 
bution from Q space, and the second summation term 
from V space of section II A . The rotation to Bogoliubov 
quasi-particles is given by the general canonical transfor- 
mation 



t 



cos ft 
e n sin ( 




(76) 



where tan2ftfc = |A|/Efc is the Bogoliubov transforma- 
tion angle. The quasi-particle annihilation and creation 
operators are indicated by d fc and 6r k . In Fig. [?| we show 
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a typical quasi-particle energy spectrum for 6 Li versus 
the single-particle kinetic energy at a magnetic field of 
B = 900 G and a temperature of T = O.QlTp. The 
figure demonstrates how well the renormalizing equa- 
tions (p2|)-(p6|) work in obtaining a cutoff independent 
energy spectrum. This is important because it implies 
that all the thermodynamics which follow will also be 
iGindependent . 

For the stationary solution the grand potential, or 
equivalently, the free energy, has indeed a minimum. 
This follows easily from setting the partial derivative 
of the grand potential with respect to <pi to zero: 
d^a/d<j>i = 0. This gives the solution 



9iV 



(77) 



which is also the stationary solution of Eq. (67). This 



equality is very useful because we can effectively elimi- 
nate the molecular field from the equations. The quasi- 
particle states are now populated according to the Fermi- 
Dirac distribution rtk = [exp(E k /ki,T) + 1] . The mean 
fields are then determined by integrating the equilibrium 
single particle density matrix elements, given by 



(27T) 



i r K 

35 J dk [(2n k - 1) cos 29 k + l] , (78) 



P 



(2^) 2 Jo 



K 



dk (2n k - l)sin26» fc . 



(79) 



Since 9 k depends on n and p, these equations require self- 
consistent solutions which are found from a numerical 
iterative method. 

In Fig. H we show a plot of the chemical potential 
as a function of temperature, for the case of 6 Li in a 
homogenous gas, at a magnetic field of B — 900 G. Fig- 
ure ^ shows the ratio of the critical temperature T c to 
the Fermi temperature Tp as a function of detuning. It 
clearly shows that there is a limiting value of T c of about 
0.5 Tf, similar to the value that has been predicted for 
40 K in Ref. ||. The BCS result for the critical temper- 
ature, given by the formula 



To 



exp 



2\a\k ? 



(80) 



has been plotted in the same graph for comparison. The 
BCS line gives a curve for T c higher than the resonance 
theory, since it does not contain the energy-dependence 
of the T-matrix. The absolute value of the scattering 
length in this magnetic field range is always larger than 
2000 ao, which implies that fef|a| > 1 — a clear indication 
that the BCS theory breaks down in this regime. 
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FIG. 8. Chemical potential as a function of temperature, 
for a magnetic field of B = 900 Gauss for 6 Li. 
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FIG. 9. Dependence of critical temperature on magnetic 
field for 6 Li, for a total density of n — 10 14 cm -3 (solid line). 
The dashed line is, for comparison, the prediction of the reg- 
ular BCS theory. 
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FIG. 10. Density profile for a gas of 6 Li atoms (solid line), 
evenly distributed among the two lowest hyperfine states. The 
temperature is T — 0.2Tf at a magnetic field of B = 900 G. 
The trap-constant is cj — 27r500 and we have a total 

number TV = 5 x 10 atoms. We compare this with a pro- 
file resulting from the same fi (but for different total number 
TV), where artificially no superfluid is present by setting the 
pairing field p equal to zero (dashed line). 

So far, our calculation has been done for a homoge- 
neous gas. We will also present results for a trapped 
lithium gas in a harmonic oscillator potential V(r) with 
a total number of TV = 5 x 10 5 atoms, similar to what 
we presented for 40 K in Ref. pi]. We treat the inhomo- 
geneity by making use of the semiclassical local-density 
approximation, which involves mainly the replacement 
of the chemical potential by a spatially-dependent ver- 
sion fi(r) = fj,— V(r). The thermodynamic equations for 
the homogeneous system are then solved at each point in 
space [p4[ . As a result, we obtain a spatially-dependent 
density distribution. At zero temperature, for a non- 
superfluid system, this gives the well-known Thomas- 
Fermi solution. For a resonance system, however, a den- 
sity bulge appears in the center of the trap, which is 
caused by a change in compressibility when a superfluid 
is present. This is shown in Figure |l^, for a spherical 
trap with a trap-constant of ui = 27r500 s _1 . This bulge 
is a signature of superfluidity and could experimentally 
be seen by fitting the density distribution in the outlier 
wings to a non-resonant system, and thus obtaining an 
excess density in the middle of the trap. For a discussion 
of the abrupt change in the compressibility see Ref. ??. 

V. FLUCTUATIONS IN THE MEAN FIELDS 
AND CROSS-OVER MODEL 

In this section we make some comments on the connec- 
tion between the resonance superfluidity theory we have 
presented and related mean-field approaches to discuss 
the cross-over of superconductivity from weak to strong 
coupling. In the mean- field theory of BEC, most often 
reflected in the literature by the Gross-Pitaevskii equa- 
tion or finite temperature derivatives, a small parame- 
ter is derived to justify the application of the theory. 
This parameter, \J na 3 may be obtained from a study 
of higher order corrections to the quasi-particle energy 
spectrum. It has been suggested that for a fermi sys- 
tem which exhibits superfluidity the small parameter is 
given by a power of kpa, and that the BCS theory breaks 
down when this parameter approaches unity. However, 
the small parameter in the theory of resonance superflu- 
idity cannot be simply a function of the scattering length 
for detunings close to resonance. This can already be 
seen from the energy-dependence of the T-matrix, which 
shows that around the Fermi-energy, the T-matrix may 
have an absolute value much smaller than at zero en- 
ergy where the scattering length is defined. Moreover, 



even right on resonance when v — and the scattering 
length passes through infinity, the T-matrix remains well 
behaved. 

Instead of calculating the small parameter of this sys- 
tem, we choose a different approach based on cross-over 
models between BCS and BEC, formulated by Nozieres 
p3[ , and later expanded upon by Randeria (34|. In 
the regular BCS theory for weakly coupled systems the 
value of the critical temperature is given by the expo- 
nential dependence in Eq. (p0|), but for strongly coupled 
systems this model results in a logarithmically diver- 
gent prediction for T c . The parameter (fcpa)" 1 is usu- 
ally taken to describe the crossover from the weak cou- 
pling Bose limit((fcFa) _1 — > — oo) to the strong coupling 
([k-pa)~ x — > +oo) BEC limit. The unphysical divergence 
in T c occurs because the process which dominates the 
transition in the weak coupling regime is the dissocia- 
tion of pairs of fermions. For a strongly coupled sys- 
tem, however, the fermions are so tightly bound that 
the wave functions of pairs of atoms begin to overlap, 
and the onset of coherence is signaled by excitations of 
the condensed state, which occurs at a temperature well 
below the dissociation temperature of the Cooper pairs. 
Thus, when moving from weak to strong coupling, the 
nature of the transition changes from a BCS to a BEC 
type mechanism. An explicit inclusion of the process of 
molecule formation, characterized by the detuning, res- 
onance width, and resonance position, will allow us to 
move from one regime to the other. 

The lowest order correction which connects between 
BCS and BEC type superconductivity can be made by 
augmenting the density equation to account for the for- 
mation of pairs of atoms. This is done by using the ther- 
modynamic number equation TV = ~^§^-, with <&g the 
total thermodynamic grand potential 

<$> G =$° G -k b TX (lAqi lnr(q,iq l ). (81) 

The term $ G is a grand potential that does not include 
the quasi-bound molecules and results from regular BCS 
theory. Retaining only this term yields a theory which 
can only account for the free and scattered fermionic 
atoms which contribute to the fermion density, therefore 
the theory breaks down if a sizeable number of bound 
states are formed. In the extreme limit of strong cou- 
pling, <£>ji. becomes negligible and equation (^Tj) just re- 
duces to the thermodynamic potential of an ideal Bose 
gas. In this regime, the theory predicts the formation 
of a condensate of molecules below the BEC transition 
temperature. 

The function r(q, iqi), which is a function of momen- 
tum q and thermal frequencies iqi, is mostly negligible for 
a weakly coupled system and has little affect on the value 
of T c in this regime. It allows for the inclusion of the low- 
est contributing order of quantum fluctuations |33|,|34| by 
means of a general inclusion of mechanisms for molecular 
pair-formation. In the resonance superfluidity model, a 
similar term is present due to the formation of bosonic 
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molecular bound states <pi , and prevents the critical tem- 
perature from diverging. When the coupling increases, 
the formation of molecules adds significantly to the total 
density equation in both the cross-over models of super- 
conductivity and in the theory we have presented here. 
Moreover, the inclusion of the molecular term allows for 
a smooth interpolation between the BCS and BEC lim- 
its. This is clearly a substantial topic in its own right, 
and will be addressed further in a future publication [|35[ . 
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VI. CONCLUSIONS 

We have shown that it is possible to derive a mean-field 
theory of resonance superfluidity, which can be applied to 
ultra-cold Fermi gases such as 6 Li and 40 K. The Hamilto- 
nian we use treats the resonant states explicitly, and au- 
tomatically builds the coupled scattering equations into 
the many-body theory. With a study of analytical scat- 
tering we have shown that these scattering equations can 
completely reproduce a full coupled channels calculation 
for the relevant energy regime. The energy-dependence 
of the s-wave phase shifts can be described by a small 
set of parameters which correspond to physical proper- 
ties, such as the non-resonant background value of the 
scattering length, and the widths and detunings of the 
Feshbach resonances. Close to resonance, we predict a 
large relative value of 0.5 Tp for the critical tempera- 
ture. The particular resonance under study for 6 Li oc- 
curs in the (1/2, 1/2) + (1/2, —1/2) collision channel, and 
has its peak at Bq = 844 G, and a width of about 
AB w 185 G @Q. This large width translates into 
a large magnetic field range where the critical tempera- 
ture is within a factor of two from its peak value. This 
range is, for comparison, much larger than for 40 K. For 
6 Li there are also two other Feshbach resonances, one in 
the (1/2,-1/2) + (3/2,-3/2) state and another in the 
(1/2, 1/2) + (3/2, -3/2) state. They result from coupling 
to the same singlet bound state, and occur at field val- 
ues of about Bq = 823 G and Bo = 705 G, and have a 
similar width to the (1/2, 1/2) + (1/2, —1/2) resonance. 
The disadvantage of these resonances, however, is that 
the atoms in these channels suffer from dipolar losses, 
which are also resonantly enhanced. Three-body inter- 
actions will be largely suppressed, as asymptotic p-wave 
collisions will give very little contribution in the tem- 
perature regime considered (an s-wave collision is always 
forbidden for at least one of the pairs). From a study 
of cross-over models between BCS and BEC we find no 
indication of break-down effects of the applied mean-field 
theory. 
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